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Abstract 

The recent suggestion that dramatic changes may occur in the Hfetime of alpha and beta decay 



, when the activity, in a pure metal host, is cooled to a few Kelvin, is examined in the light of 



published low temperature nuclear orientation (LTNO) experiments on such sources cooled to 
as low as 25 mK, with emphasis here on alpha decay. In LTNO observations are made of the 



o 

D . anisotropy of radioactive emissions with respect to an axis of orientation. Correction of data for 



O I decay of activities in metallic samples held at temperatures at and below 1 Kelvin for periods of 

^ ' several days has been a routine element of LTNO experiments for many years. No evidence for any 



change of half-life on cooling, with an upper level of order 1%, has been found, in striking contrast 
to the predicted changes, for alpha decay, of several orders of magnitude. The proposal that such 
dramatic changes might alleviate problems of disposal of long-lived radioactive waste is shown to 
be unrealistic. 
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I. INTRODUCTION 



In a recent paper ideas were presented which led to the suggestion that the half-hves 
of alpha emitters, particularly those of lower decay energy, would be strikingly reduced if 
the activity was in a metallic environment and cooled to temperatures of a few Kelvin as 
compared to 'normal' decay in insulating materials at room temperature The authors 
surmised that, through the influence of the conduction, 'valence' or 'free', electron potential 
near the nucleus, the alpha decay barrier would be effectively lowered and thus cooling of 
long-lived radioactive waste might lead to their much more rapid decay and safe disposal. 
No experimental evidence was offered in support of these suggestions. It is the objective 
of this paper to report aspects of already published experiments of significance to these 
suggestions, clearly indicating that the predictions are far from experimental reality. The 
experiments concern the decay of relatively long lived ^^^Rn, ^^^Ra and ^^^Ac and their 
daughter decay chains implanted in an iron metal foil and studied at temperatures down 
to below 20 mK. The original objectives of this work were to study the polarization of the 
decaying isotopes through measurement of the angular properties of their alpha and gamma 
emissions - the method of low temperature nuclear orientation (LTNO) - leading to results 
on the hyperfine interactions of the isotopes involved, their magnetic dipole moments and 
the magnetic hyperfine fields they experienced as implants in the ferromagnetic iron lattice 
aaa. However the .e.u.ts throw direct hght upon the sugge.t,o,rs ,„ade by Kettner et al 
in Ref. p|. Brief mention of evidence from experiments involving beta decay is made in the 
discussion. 

This paper presents experimental evidence only, without detailed comment on the theo- 
retical considerations on which the proposals of Ref. |l| were based. We note however that 
(i) serious objections to the underlying ideas on the effect of screening on alpha particle 
Coulomb barrier penetration have been raised elsewhere |5.] and (ii) that the suggestion 
that cooling to low temperatures might be worthwhile depends upon treating conduction 
electrons as classical gas with mean energy proportional to temperature T. Furthermore, 
the idea that screening can take place by electrons at a Debye radius much smaller than 
the K-shell mean radius , falling towards nuclear radii at and below IK, contravenes 
basic quantum mechanics. There are simply no states available to construct such confined 
electron wavef unctions. If there were they would be occupied by the atomic electrons. More 
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comments are given in the discussion. 

The experiments which provide the evidence given here were performed at ISOLDE, 
CERN and Oxford over a period of several years. Other LTNO groups, in particular those 
from Leuven and Bonn, have published work on implanted alpha emitters cooled to mil- 
likelvin temperatures in metals, containing similar evidence (see e.g. 0, ^,[9] and references 
therein) . 

The basic ideas of the proposed enhancement mechanism are given in Sec. HIl which is 
followed in Sec. IIIII by a survey of the ways in which low temperature nuclear orientation 
experiments are sensitive to source lifetime and decay correction. Sec. II VI gives details of the 
preparation of the implanted metallic samples used in the work referred to and Sec. |V] in- 
troduces model calculations based on three possible scenarios regarding the correct lifetimes 
to use in analysis of the experiments. Sec. |Vl] gives an extended analysis of predictions and 
experimental results for the three decay chains headed by ^^^Rn, ^^^Ra and ^^''Ac. Discussion 
and conclusions are presented in Sec. IVIII 

II. BASIC IDEAS OF THE PROPOSED ALPHA DECAY ENHANCEMENT 
MECHANISM 

From Ref. [l| the enhancement of alpha decay rates in a metal is associated with the 
presence around the metallic ion of a screening charge, due to free electrons, which is taken 
to reduce the barrier height by the screening energy Ud- According to this picture, the 
reduced barrier height enhances the rate of alpha tunneling through the barrier to separate 
from its parent nucleus, thus reducing the lifetime. Since ionic charge screening by free 
electrons occurs only in metals, it is predicted that large differences of alpha decay lifetime 
will be observed depending upon whether the decay takes place in an insulator or a metal. An 
expression is given for Ud in terms of the charges of the alpha particle Zq, and the daughter 
nucleus Zt, the temperature of the lattice electrons and a parameter Ue(d+d) which relates 
to a d + d fusion reaction in a metal 

/290\ 

U^ = Z^ZMd + d)2go[-jr) (1) 

where T is the absolute temperature of the metal in which rtie decay takes place. A 'typical' 
value of Ue(d+d)29o is given to be 300 eV at 290 Kelvin [1]. The temperature dependence 
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arises from taking the Drude model for the 'quasi-free' valence electrons as having an (av- 
erage) kinetic energy of 0.5 kT. Based on this expression it is suggested that by cooling the 
activity in a metal sample to liquid helium temperature (4.2 K), the half-life of ^^°Po can be 
reduced from 138 d to 0.5 d and of ^^^Ra from 1600 y to 1.3 y, factors of order 1000, in each 
case through a reduction of the effective barrier potential by about 420 keV (see Ref. 

The phenomenon of ionic charge screening in metals is well established and its conse- 
quence is giving rise to small energy shifts, of order several hundred electronvolts, for cross 
sections in low energy reaction studies in metallic targets, is fully documented in Ref. l(] 
and other sources. The issue in this paper is to explore whether, through a strong depen- 
dence of the effect on temperature, charge screening can lead to dramatic consequences at 
temperatures of order a few degrees and below. Three scenarios are considered and their 
consequences compared with existing experimental evidence. They are: 

• The Full effect: the suggested mechanism operates fully at all temperatures. 

• The Partial effect: an attenuated version of the mechanism occurs, leading to a 5% 
reduction in all alpha decay half-lives at 1 K, which however obeys the predicted Ud 
temperature dependence. 

• The Null effect: all effects of screenimg on half-life are negligibly small so that 'normal' 
half-lives may be used. 



III. LOW TEMPERATURE NUCLEAR ORIENTATION EXPERIMENTS AND 
THEIR SENSITIVITY TO HALF-LIVES. 

Low temperature nuclear orientation experiments involve cooling radioactive isotopes in 
samples in which they experience strong hyperfine interactions to temperatures at which 
they become, in the case of a magnetic interaction, polarised along the direction of the 
magnetic field. The appropriate temperature range is typically between 1 and 100 mK for 
activities in, for example, a ferromagnetic iron foil sample. For relatively long-lived activities 
such as are concerned here, the experiments are performed with the sample soldered to the 
cold finger of a dilution refrigerator and involve initial cooling, by stages, to liquid nitrogen 
(77 K) and liquid helium (4.2 K) temperatures, each for about 24 hours to allow for cooling 
and establishment of good vacuum conditions (less then 10~^ mbar), before reaching 1 
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K attached to a pumped liquid helium reservoir. At 1 K temperature the nuclei remain 
unoriented and 'warm' counts, to be used for normalisation of the measured anisotropies at 
lower temperatures, are taken for several hours. Circulation of the ^He/^He mixture in the 
refrigerator is started and after a few hours the sample temperature enters the millikelvin 
range. It can be kept there for many days, but often circulation is stopped to allow return 
to 1 K to take more 'warm' counts in order that the decay of the sample can be fully 
allowed for in analysis. In the orientation measurements themselves the basic observable is 
the anisotropy of angular distribution of emission of particles or radioation from an oriented 
source 

A={W{e,T)-l)% (2) 

where W{9,T) is the angular distribution function at angle 9 with respect to the axis of 
orientation and temperature T. It is derived form experimental count rates by 

Ar(^) and A^(™) are decay and dead-time corrected count rates in a detector when the source 
in 'cold' (oriented at millikelvin temperature) and 'warm' (unoriented at ~ 1 K or above), 
respectively. The decay correction, fundamental to these measurements, is directly relevant 
to the present work, as it involves the decay constant of a isotope in question through a 
correction factor 

A5exp(+At) 



1 — exp(— A5) 

where A=ln(2)/Ti/2 is the decay constant, 5 is the real time of the counting period and t is 
the difference between the arbitrary normalisation time and the start of the measurement. 
After this correction is applied, each count rate can be compared to rates in spectra taken 
at different times as though they were taken at the same instant in time. Physical param- 
eters are extracted from measured anisotropies by comparison of the experimental angular 
distribution and the theoretical expression 

^max 

W{e, T) = 1 + / £ 5k(T)[/kAkQkPk(cose) (5) 

k=2 

where Bk, Uk and Ak are the usual orientation, deorientation and angular distribution co- 
efficients, Qk are the detector solid angle corrections, Pk are Legendre polynomial and / is 
the fraction in good s.es y . Depending on .Ke na..e of a p„.tien>. expen.en, .elected 
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information contained in one of the coefficients can be obtained by such a comparison, 
assuming all the other coefficients are known. In the work referred to in this paper, the 
orientation coefficients B^, dependent on the product of magnetic dipole moment and the 
hyperfine field acting on the studied nuclei in the sample, were used to extract the mag- 
netic dipole moment of oriented states and compare them with either the same moments 
obtained by a different technique or with systematics and theory. It is obvious that should 
the experimental anisotropy be determined incorrectly, i.e. if the correction for the radioac- 
tive decay, made by using decay constants known at room temperature, would not apply at 
temperatures of 1 K and lower, the extracted magnetic dipole moment would be wrong. 

Additional sensitivity to the appropriateness of the decay correction arises from other 
detailed properties of the measured anisotropic distributions. The first sensitive feature is 
that, at the lowest temperatures, the anisotropy must become independent of temperature 
as the nuclei become fully polarized. Since there is a time sequence in such measurements 
as the sample temperature is changed, any incorrectness in the source decay correction will 
impose a false slope or scatter upon the extracted temperature dependence. The second 
sensitive feature is that the anisotropy involves an increase in observed intensity in some 
directions to the orientation axis, and in others a decrease, the actual quantitative changes 
being determined by well established theory of the distribution. The application of an 
incorrect adjustment for source decay will falsify the relationship between observations in 
detectors in different directions to the axis. Thus, if the decay correction assumes too short 
a lifetime, the 'warm' unoriented distribution will be taken to fall too fast, so that any 
observed increase in intensity, measured at a time later than the 'warm' distribution, will be 
analysed as too large relative to its true value by a certain factor, and any simultaneously 
observed decrease would be analysed as smaller than its true value, by the same factor. 
These failures would lead to distortion of the relation between the analysed anisotropics in 
the two directions, which can usually be readily identified. 

A separate aspect of the LTNO technique with relevance to hfetime is the fact that for the 
observed anisotropics to be characteristic of the hyperfine interaction in a specific nuclear 
state, usually a ground state, the lifetime of that state must be long enough for the nuclei 
to reach thermal equilibrium with the metallic lattice and achieve a degree of polarisation 
characteristic of that temperature and interaction strength. The time for this to be reached, 
the nuclear spin-lattice relaxation time, Ti, has been established in implanted, diffused and 
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co-melted metallic samples as following the Korringa conduction electron mechanism. In 
iron samples of the type discussed in this paper Ti's have been measured to range from 
many hours to below 100 ms at millikelvin temperatures. The relevant parameter is the 
ratio Ti/ Ti/2- There is abundant evidence that were a lifetime to become shorter than 0.1 
ms thermal equilibrium will not be achieved. In situations involving a decay cascade, such 
as are described in this paper, a degree of polarisation may be 'inherited' from a precursor 
state higher in the decay chain. Reference to situations of this kind are made in the text 
where relevant, as the degree of equilibration will be affected if lifetimes were to be severely 
reduced on cooling. 

Finally, as a matter of course, OLNO experiments include measurement of the unoriented 
intensity at both the start and the end of the experiment, and quite often also at an inter- 
mediate time, when the temperature is at close to 1 K. It is a standard requirement that 
these normalisation measurements must, after due correction for decay of the source, yield 
a consistent value for the reference 'warm' intensity to be used in evaluation of all measured 
anisotropics. If the normalisation should not be consistent this would be immediately ap- 
parent and present a challenge to the experimenters. Thus we have, in effect, a measure of 
the source decay at 1 K over the period of the experiment, which is typically at least several 
days. 

Each experiment has individual sensitivity to sample half-life and the cases of the three 
alpha decay chains described here are discussed in detail in Sec. IVD 

IV. IMPLANTED METALLIC SAMPLE PREPARATION FOR LTNO ALPHA 
ANISOTROPY EXPERIMENTS. 

For the experiments described here, alpha active samples of ^^^Ra, ^^^Ra and ^^^Ac were 
prepared at the ISOLDE isotope separator facility, CERN, the parent activities being im- 
planted at 60 keV into prepared 99.999% pure iron foils which have been carefully surface 
polished. The implant dose was ~2-10^'^ for mass A= 225 and A=227 and ~2-10^^ for 
mass A=224. Experience of implantation of a wide range of elements has shown that the 
implantation energy, 60 keV, is amply sufficient to allow the implanted ions to penetrate 
the omnipresent thin oxide layer on the foil surface and come to rest in undisturbed pure 



iron. Estimates based on calculated projected ranges llj] show that an ion with a nucleus of 
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mass A~200 accelerated through 60 keV will penetrate ~14 nanometers, or about 75 atomic 
layers in iron. Allowing for an oxide layer of approximately 5 nanometers the average 
depth of the initially implanted alpha emitters is ~45 atomic layers of pure iron. In samples 
where there is an alpha decay cascade, recoil following alpha emission alters the initial depth 
distribution. The ~100 keV nuclear recoil energy can cause ejection of the daughter from 
the sample, but, on average, the mean implant depth of those nuclei remaining in the pure 
iron, about 75% of the initial quantity after a three-stage cascade, is increased considerably 
compared with the initial depth [l^. Many experiments have demonstrated that this im- 
plantation energy is adequate to describe the environment of the implants as fully typical of 
the metal. In the context of low temperature nuclear orientation evidence comes from the 
fact that Nuclear Magnetic Resonance (NMR) of the implants has been observed at the same 
frequencies in implanted samples as in samples prepared by thermal diffusion to depths of 
~1 micron [l^. Additional evidence of the implants being in true metallic lattice is the fact 
that the implanted nuclei have been observed to relax in on-line implantation experiments 
with time constants Ti described fully by the Korringa conduction electron mechanism jisl . 

V. PREDICTIONS OF ALPHA ACTIVITY CHANGES ON COOLING 

As mentioned is Sec. [Tll we consider detailed estimates of effects of proposed lifetime 
reductions in three possible scenarios to be referred to as the Full, Partial, and Null cases. 
In all cases the lifetime reduction estimates, based on the expression for the potential Ud, as 
given in Ref. jl| and quoted above, can be readily estimated, without detailed calculation, 
3y__taking the systematic half-lives of allowed alpha emitters versus decay energy (see e.g. 



16|) to obtain d(log Ti/2)/dEQ, as a function of the alpha particle decay energy, 



TABLE I: Predicted alpha decay lifetime reductions and alpha fraction changes 
based upon Ref. [l| 



Parent Normal alpha Alpha Normal 


Ud 


HALF-LIFE 


PREDICTED ALPHA 


PREDICTED 


Nucleus decay half-life Energy fraction 


(keV) 


reduction 


FRACTION 


DECAY HALFLIFE 


AT 290 K * (keV) (%) AT 290 K 




FACTOR AT IK 


(%) AT IK ** 


{%) AT IK *** 


A=224 decay chain, headed by alpha emitter ■^■^'*Ra 










224Ra 3.66 d 5686 94 


898 


6.3 X 10* 




4.7 S 


5449 5.5 


898 


1.3 X IQS 
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TABLE I: (continued) 



Parent 


Normal alpha 


Alpha 


Normal 


Ud 


HALFLIFE PREDICTED ALPHA 


PREDICTED 


Nucleus 


DECAY HALFLIFE 


Energy fraction 


(keV) 


REDUCTION 


FRACTION 


DECAY HALFLIFE 




Ai 21)0 K ' 


(Ki;\') 


('a) ai 2!)() K 




FA( KJH AI IK (\ 


/;) A'l IK 


{'A) A'i iK 


220Rn 


55.6 s 


6228 


100 


878 


1.2 X 10* 




4.7 ms 


216po 


0.15 s 


6779 


100 


857 


2.7 X 10^ 




55 X 10-'^ s 


2123; 


60.6 m 


6090 


27.2 


847 


1.2 X 10* 


26 


0.82 s 






6051 


69.9 


847 


1.3 X 10* 


74 




predicted 


missing transitions if enhanced alpha decay of weakens beta decay branch 




212Po 


0.30 X 10-s 


8784 


100 


not seen if beta decay of 2123; ceases 




A=225 decay chain, headed by ^^^Ra. >99% beta <1 x 10 


-*% alpha Ti/2 14.i 


i d 




225 Ac 


10.0 d 


5830 


50.7 


889 


3.8 X 10* 


46 


20.6 s 






5794 


24.3 


889 


4.3 X 10* 


25 








5732 


10.1 


889 


5.0 X 10* 


12 








5637 


4.4 


889 


6.5 X 10* 


7 




22lpr 


4.8 m 


6341 


83.4 


867 


7.7 X 10^ 


77 


34 ms 






6126 


15.1 


867 


1.3 X 10* 


23 




2i7At 


32.3 ms 


7067 


99.9 


867 


1.8 X 103 


100 


18 X 10-6 s 


213Po 


4 X lO"'' s 


8375 


99.9 


857 


3.3 X 102 


100 


1.2 X 10"* s 


predicted additional transitions via 


enhanced alpha 


decay 


of 225Ra 






22BRa 


>4 X lO^y 


4950 


not seen 


898 


5.0 X 10^ 


30% of all 


28 d 


221 Rn 


125 m 


6035 


strongest 


878 


2.0 X 10* 


strongest 


0.38 s 


217po 


<10s 


6537 


>95 


857 


4.4 X 103 


100 


< 2 ms 


predicted additional transitions via 


enhanced alpha 


decay 


of 213Bi 






213Bi 


38 h 


5869 


93 


847 


2.1 X 10* 


85 


5.9s 






5549 


7 


847 


4.8 X 10* 


15 




A=227 decay chain headed by 227 A 


c , normal decay; 98.6% beta, 1.4% alpha, 


Ti/2 21.8 y 




227Th 


18.72 d 


6038 


24.5 


918 


3.2 X 10* 


19.7 


40.7 s 




100% alpha 


5979 


23.4 


918 


3.5 X 10* 


20.5 








5756 


20.3 


918 


6.8 X 10* 


34.7 








5713 


4.9 


918 


7.6 X 10* 


9.3 








5709 


8.2 


918 


7.6 X 10* 


15.6 
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TABLE I: (continued) 



Parent 


Normal alpha 


Alpha 


Normal 


Ud 


HALFLIFE 


PREDICTED ALPHA 


PREDICTED 


Nucleus 


DECAY HALFLIFE 


Energy 


FRACTION 


(keV) 


REDUCTION 


FRACTION 


DECAY HALFLIFE 




AT 290 K * 


(keV) 


(%) AT 290 K 




FACTOR AT IK 


(%) AT IK ** 


(%) AT IK *** 




11.44 d 


5748 


9.1 


898 


5.4 X lO* 


7 


14.0 s 




100% alpha 


5717 


53.7 


898 


6.1 X 10* 


46.7 








5607 


26 


898 


7.9 X 10* 


29.1 








5539 


9.1 


898 


9.9 X 10* 


12.7 








5434 


2.3 


898 


13.5 X 10* 


4.4 




2l9Rn 


4 s 


6819 


81 


878 


3.1 X 103 


62.9 


1.1 ms 




100% alpha 


6553 


11.5 


878 


5.3 X 10^ 


16.6 








6425 


7.5 


878 


7.2 X 10^ 


14.7 




215Po 


1.8 ms 


7386 


100 


857 


1.0 X 10^ 


100 


1.8 X lO"'' s 




2.15 m 


6623 


84 


847 


3.3 X 10^ 


71 


33 ms 






6279 


16 


847 


7.1 X 10^ 


29 




predicted additional transitions via enhanced alpha 


decay 


of 227 Ac 






Ac 


21.77 y 


4950 


47 


889 


5.7 X 10" 


45.4 


1.00 d 






4938 


40 


889 


5.7 X 10^ 


38.8 








4870 


6.1 


889 


7.5 X 10^ 


7.8 








4853 


3.7 


889 


7.5 X 105 


4.7 




223pj, 


21.8 m 


5340 


strongest 


867 


1.3 X 10^5 


100 


3.4 m 


219 At 


0.9 m 


6275 


strongest 


867 


4.0 X 103 


100 


13.5 ms 



* allowing for beta decay branching 

** after renormalisation of alpha branching 

***neglecting any change in beta half-life 



This simple approach is vahd as the half-hfe adjustment relates to the barrier penetra- 
tion factor, not to the pre-formation of the alpha particle within the parent nucleus. This 
method, with estimates of direct from Eqn. ([1]) and the typical room temperature value 
of Ue(d+d)29o given in Ref. has been used to obtain estimated half-life reductions, see 
Table [U for the principal alpha decay transitions involved in the three isotope chains rele- 
vant to this paper at 1 K. These are examples of the Full scenario. Estimates of changes at 
different temperatures (77 K, 4.2 K and 25 mK) for the Full scenario have been made by 
the same method. For the Partial scenario it is assumed that the lifetime change is -5% at 
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FIG. 1: ^^^Ra decay chain; f3 1/a branching ratio from ^^^Bi predicted to weaken. 

1 K and the simple theory have been used to give the value of Ud to which this reduced 
effect would correspond. It has been then further assumed that, to estimate the effects at 
different temperatures, the temperature dependence follows the predicted T^^/^ relationship 
of Eqn. [TJ The Null scenario makes no allowance for any half-life change with temperature. 
This is the expectation were the screening electrons to obey Fermi-Dirac statistics and has 
been used implicitly in all published OLNO experiments. 

VI. EXPERIMENTAL BEHAVIOUR OF COOLED ALPHA SOURCES IM- 
PLANTED INTO IRON, COMPARED WITH THE FULL, PARTIAL AND NULL 
SCENARIOS 

The decay chains of the three implanted isotopes listed above are shown in Figs. [T][5] and 
E] and details of the lifetimes and stronger alpha decay transitions in each decay are given 
in Table [H All half-lives are taken from [tt] . The consequences of the predicted half-life 
changes are qualitatively different for the three chains as outlined below. 
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50 100 150 200 250 

Time [hours] 



FIG. 2: Relative activity for selected members of the Ra decay chain in the Full and Null scenar- 
ios (top panel) and sample temperature (bottom panel) as a function of time after implantation. 

A. The ^^"^Ra decay chain. 

The implanted isotope, ^^^Ra is an alpha emitter (3.66 d) and a single sample was studied 
during an experiment at the NICOLE dilution refrigerator at ISOLDE, CERN, which lasted 
for about 8 days starting one day after the sample was implanted. The time and temperature 
sequence of this experiment is shown in Fig. [2] which also illustrates the predictions based 
on the Full and Null scenarios. The Null scenario is simple in that, after the first 20 hours 
or so, all activities are predicted to follow the 3.66 d half-life of the long-lived parent. The 
Full scenario, at the other extreme, predicts that the parent activity, which starts out as the 
upper line in Fig. [21 should show two sharp step increases as the sample is cooled, first to 
liquid nitrogen temperature, 77 K, then to 4.2 K. At this temperature the predicted half-life 
is already shorter than the time which elapsed before being cooled further to 1 K, at which 
the lifetime is predicted to be 4.7 s (Table [T]) so the ^^^Ra activity falls away very steeply. 
The three prominent alpha transitions at 5686 keV (^^^Ra), 6288 keV (^^°Rn) and 6779 keV 
pi6pQ^ should thus essentially vanish before the OLNO experiment started, leaving only the 
(temporarily increased) ^^^Pb beta activity, and its daughters, present in the source, with 
alpha transitions at 6090 keV and 6051 keV, which decay with the normal ^^^Pb half-life. 
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50 100 150 200 250 300 

Time [hours] 



FIG. 3: As Fig. [2] but the Partial and Null scenarios. 

10.6 h. 

The LTNO experiment has two periods at very low temperature, taken to be for simplicity 
a steady 25 mK, with periods at IK at the start, after about 110 hours and at the end, about 
190 hours after the sample was made. It is clear that, in the Full scenario, the ^^^Bi activity 
should fall steadily throughout the experiment with the a 10.6 h half-life. 

The Partial scenario is illustrated in Fig. [3l which compares the predictions of this case 
with the Null scenario. For the Partial case there is little change until the sample reaches 
1 K at which there is a small step upward in all decay rates. A much larger step occurs on 
cooling to 25 mK, immediately in the case of ^^^Ra and delayed through the (unchanged) 
10.6 h beta decay half-life of ^^^Pb for the activities of ^^^Pb and ^^^Bi. The alpha decay 
fraction of ^^^Bi is predicted to increase from its normal value of 36% to 37.2% at 1 K and 
43.5% at 25 mK. It is clear from Fig. [3] that a set of 'warm' counts taken during the three 
periods of the experiment with the sample at 1 K would not follow the predictions of the Null 
scenario were the Partial scenario to be correct. Furthermore large distortions of the ratios 
of 25 mK activity to 1 K activity in detectors placed at angles for which the true anisotropy 
would lead to increase/decrease of counts rates would be produced by the common increases 
predicted on cooling the sample to 25 mK. 
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10 20 30 40 50 60 70 
1/T [Kelvin ^] 

FIG. 4: ^^^Ra decay chain: ^^^Bi 727 keV gamma transition anisotropy vs inverse temperature Q]- 
1. Experimental findings for the '^'^'^Ra decay chain. 

In the published experiment on this decay chain only data the 727 keV gamma transition 
from the daughter of ^^^Pb, ^^^Bi, was analysed. Fig. H] shows the experimental results 
which gave anisotropies up to about 7% at the lowest temperatures. The raw data were 
taken over a period of about 5 days and all decay corrections were made using the simple 
3.66 d half-life of ^^^Ra. The observed temperature dependence was fitted with the magnetic 
dipole moment of ^^^Bi as the unknown quantity, yielding the value | | = 0.41(5) n.m. j^, 
which is in very reasonable agreement with a subsequent laser-based measurement of 0.32(4) 



n.m. 



y. 



These observations are in complete contradiction to the Full scenario prediction 
that the the ^^^Bi activity would follow a 10.6 h lifetime throughout the low temperature 
experiment. 

As Fig. [3] shows, on even the Partial scenario, the three sets of 1 K 'warm' counts would 
have been seriously discrepant, the first and third lying close to the time dependence pre- 
dicted by the Null scenario, but the second being above them by 24%, which was not found, 
all three sets being in good statistical agreement. A second, equally unobserved, prediction 
is that the counting rates in all detectors would have been shifted by 36% on cooling to 
millikelvin temperature. We conclude that the published data provide evidence that any 
half-life change on cooling to IK is well below 5% (the Partial scenario assumption), less 
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FIG. 5: Ra decay chain with additional transitions if alpha decay enhancement were strong in 
225Ra and 2i3Bi. 

than 1 part in 10^ of the Full scenario prediction, and that a realistic upper limit of any 
effect is of order 1%. The same upper limit of 1% applies to any changes between 1 K and 
25 mK. 

B. The 22^Ra decay chain. 

The case of the ^^^Ra decay chain is basically different since the long-lived parent, ^^^Ra 
normally decays by beta decay with a half-life of 14.8 d to ^^^Ac with no alpha branch 
and is not predicted to show (major) half-life change. ^^^Ac, with 'normal' half-life 10.0 d 
then initiates an alpha decay sequence to ^^^Bi. Discussion of possible changes in the alpha 
spectrum caused by changing j3~/a branching ratios in this decay chain (e.g. ^^^Ra and 
2^^Bi) is omitted as such effects should be much stronger in the ^^'^Ac decay (see below). 

The predicted behaviour on cooling under the three scenarios is illustrated in Fig. [61 As 
shown in the figure, this sample was the subject of three experiments. The first two took 
place at ISOLDE, CERN, using the NICOLE dilution refrigerator. Of these the first started 
about 24 hours after the sample was made and continued for about 350 hours whilst the 
second started after about 750 hours and ended after about 820 hours. The sample was 
then taken to Oxford and the last experiment was run with the sample between about 1300 
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FIG. 6: Relative activity for ^^'^Bi and ^^^Ra in the ^^^Ra decay chain in the Full, Partial and Null 
scenarios (top panel), and sample temperature (bottom panel) as a function time after sample 
implantation. 

and 1500 hours old. For each experiment the temperature sequence was similar to the ^^^Ra 
chain experiment described above. In Fig. [6] the decay of the parent activity is a straight 
line with 'normal' half-life of 14.8 d, whilst the 'normal' (10.0 d) decay of ^^^Ac and its much 
shorter lived daughters should follow the rising and then falling smooth curve, which is the 
Null scenario prediction for this decay chain. The Full scenario prediction is that, in each 
of the three experiments, the ^^^Ac activity, and that of its immediate daughters, should 
show a large increase to well above the 'normal' growth curve at 77 K, which disappears 
over about 30 hours when the sample is held at this temperature. Before equilibrium is 
re-established however, the temperature has been reduced to 4.2 K, leading to a second 
increase in count rate, which died away within less than hour. The activity then settles 
into simple equilibrium with the decay of ^^^Ra parent, with a 14.8 d half-life. Subsequent 
cooling to 1 K and to millikelvin temperature, and return to 1 K during the experiment, has 
little effect in this scenario since there is very little activity 'held-up' in the decay chain for 
such short daughter half-lives. 

The ^^^Bi beta activity will show similar strong short term increases, which recover with 
close to its 'normal' half-life of 45.6 m before also coming to equilibrium with the 14.8 d 
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decay, showing no effect of the intervening activities. It is predicted to decay with a 14.8 d 
half-hfe for the duration of the sub-4.2 K periods of all experiments. 

After each experiment at low temperature, when the sample is returned to room temper- 
ature, the activities are predicted to fall dramatically then recover towards the Null scenario 
curve. This recovery has the time constant of the 'normal' 10.0 d ^^^Ac half-life and so should 
take about three weeks, longer than the interval between the first and second experiments. 

In Fig. [6] the Partial scenario is seen as producing relatively minor, but still distinctly 
visible, changes compared to the Null scenario. Any effect before reaching 1 K is small, but 
the Partial scenario gives a predicted activity increase of 28% in all detectors on cooling 
between 1 K and 25 mK. 

1. Experimental findings for the '^'^^Ra decay chain. 

In the first two experiments on this chain only gamma decay data were taken, whilst in 
the third alpha decay spectra were also observed down to millikelvin temperatures. The 
published record shows anisotropies measured on the alpha transitions from ^^^Ac, on both 
alpha and gamma transitions from ^^^Fr and from the alpha decay of ^^^At, all precursors 
of the decay of ^^^Bi Data for the 292 keV gamma transition in the decay of ^^^Bi is 
shown in Fig. [71 which gives results for detectors along (0°) and normal to (90°) to the axis 
of nuclear orientation. The data throughout the decay sequence were treated with decay 
corrections based on the 'normal' lifetimes of the longer lived components of the chain, ^^^Ra 
(14.8 d) and ^^^Ac (10.0 d), following the Null prediction. Several features of all three 
experiments give clear evidence contrary to the Full scenario. The most obvious is that the 
time dependence of the 1 K 'warm' counts was in agreement, within small ~1% statistical 
variation, with the Null prediction, in all experiments even though the time variation during 
each experiment differed depending upon the age of the sample (see Fig. [6]). Secondly, very 
different temperature dependence was shown for the alpha anisotopies in the decays of ^^^Fr 
('normal' Ti/2= 4.8 m) and ^^''At ('normal' Ti/2= 32 ms). This is only possible, as outlined 
in Sec. IIIH if there is relaxation of these nuclei to achieve thermal equilibrium with the 
surrounding metallic iron lattice within their lifetimes. In the Full scenario, from Table. [H 
the lifetimes of ^^^Fr and ^^^At become 34 ms and 18 fis at 1 K. Both become less than 1 
fis at 25 mK, far shorter than the relatively well known Korringa mechanism would be able 
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FIG. 7: ^^^Ra decay chain: ^^^Bi 292 keV gamma transition axial (0°) and equatorial (90°) 



anisotropy vs inverse temperature 
to relax to equilibrium 



15j . Thus these activities should not become appreciably oriented 
unless fed from an already oriented source. Thus ^^^Fr and ^^^At should both show inherited 
from the parent ^^^Ac. However, ^^^Ac has spin 1/2 and for spin 1/2 isotopes no anisotropy 
is observed in emissions which conserve parity. Thus no orientation detected by alpha or 
gamma decay should be observed in emissions from either 22ipr or ^^^At in the Full scenario 
in direct contrast with experiment Sj. 

To place an upper limit on any possible effect, predictions of the Null and Partial scenarios 
are again compared (Fig. [7]). The most important feature of the Partial scenario is that there 
is a clear discontinuity in activity as the sample is cooled to 25 mK from IK which, as detailed 
above, must distort the 0^/90'' anisotropy relationship. The data on the 292 kev transition 
are of high statistical accuracy (~0.5%) and show excellent agreement with theory, having 
the correct ratio to one another and both tending to a constant value as the temperature 
approaches 16 mK (1/T = 60 K^^). Any deviation greater that 2% would be revealed, as 
compared with the predicted shifts (which are both upward, i.e to increase the 90° effect 
and reduce that at 0°) of ~25% in the Partial scenario. 

The observed anistropies show no sign of such distortions. Both alpha and gamma results 
were analysed to give the magnetic dipole moment of ^^^Bi. The value obtained, \ fi \ = 
3.89(9) n.m. [S], is in very good agreement with a subsequent laser-based measurement of 
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FIG. 8: ^^''Ac decay chain with additional transitions if alpha decay enhancement were strong in 

227Ac. 

3.716(7) n.m. Q. 

On this basis the evidence from the ^^^Ra mass chain is that there is no half-hfe effect 
above 0.5%. 



C. The 'Ac decay chain. 

The 22^Ac chain is shown in Fig. [71 including predicted changes associated with the Full 
scenario described below. The 'normal' sequence is by dominant (98.6%) beta decay to 
22^Th, with half-life 21.8 y, followed by the strong alpha sequence ^^^Th - ^^^Ra - ^^^Rn - 
2i5pQ _ 2iip|^^ There is a weak (1.4%) alpha decay branch to ^^'^Fr, which returns to the main 
sequence through (99%) beta decay of ^^^Fr to ^^^Ra. A second weak (0.005%) alpha branch 
{Ea ~ 5300 keV) exists from ^^^Fr to ^^^At which then decays further by dominant alpha 
decay to ^^^Bi. The Null scenario prediction is that this decay follows a time dependence 
basically determined by the time for the two longer-lived alpha decays, ^^''Th (18.7 d) and 
22^Ra (11.4 d), to come into temporal equilibrium with the 21.8 y ^^''Ac activity. Fig. [9] 
shows that this was complete before the experiment, which took place approximately 150 
days after the sample was implanted. Both alpha and gamma transitions from all elements 
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of the decay chain were observed and alpha spectra are shown at room temperature (Fig. [T0|) 
and at 1 K and ~10 mK (Fig. [TTl). 

The Full scenario estimated half-life changes from Table [T] alter this sequence fundamen- 
tally at 1 K ( and there will be less striking changes at higher temperatures). Through alpha 
decay enhancement, the I3~ /a branching ratio for ^^'^Ac is altered from 70 to 1.5 x 10~^ so 
a decay to ^^^Fr should become dominant. Furthermore the subsequent decay of ^^^Fr also 
becomes dominated by alpha decay (86%) with beta decay reduced to 14%. Thus, as the 
sample is cooled, not only should the total activity show marked increases but very strong 
feeding of ^^^At and ^^^Bi should be observed before the source loses activity with a half 
life of about 1 day, rather than the normal 21.8 y. A further major prediction of the Full 
scenario for this decay chain is that the alpha decay fine structure observed in the decays 
of ^^^Th and ^^^Rn will be seriously altered, the higher energy alpha transitions becoming 
weaker relative to the lower energy transitions. The calculations presented in Tabled predict 
a doubling at 1 K of the intensity of the 5756 keV transition in the decay of ^^^Th relative to 
the 6038 keV transition. In ^^^Ra the relative intensity of the 5434 keV alpha branch should 
increase from 2.3% to 4.4% whilst the 5717 keV branch is predicted to weaken from 54% to 
47%. These combined changes would seriously modify the appearance of the spectral region 
between 5300 and 6100 keV in the spectrum from this sample at 1 K as compared to room 
temperature. Cooling to millikelvin temperatures would, in this scenario, lead to further 
major changes in half-life and fine structure relative intensity. 

The Partial scenario predictions for the experiment performed on the ^^''Ac decay chain 
are shown in FigO The small change in the ^^^Ac (3~ /a branching ratio has been neglected 
so that the experiment starts close to the temporal equilibrium situation of the Null scenario. 
Cooling to 4.2 K produces ~1.6% increase in all activities, with a further 1.2% increase at 
1 K and a much larger (~30%) increase on cooling to 25 mK. These increased activities die 
away only slowly during the 1 K and 25 mK sections of the experiment, with the exception 
of ^^^Th which shows a decay with the reduced lifetime of ~13 d throughout in this period. 
This decay is far faster than the Null scenario and has direct consequence for the sequence 
of 1 K 'warm' counts on alpha transitions from this isotope. 



20 



50 



45 



a 40 



OS 



Ac Null 

7 

Th Partial 

3 

Ra Partial 
'Bi Partial 



35 



V 




3650 



3700 



3750 3800 
Time [hours] 



3850 



3900 



FIG. 9: Relative activity for selected members of the ^^''Ac decay chain in the Null and Partial 
scenarios (top panel), and sample temperature (bottom panel) as a function of time after implan- 
tion. 

1. Experimental findings for the '^'^'^ Ac decay chain 



In the data analysis, all decay corrections were made with the 21.8 y ^^'^Ac half-life, the 
Null scenario, and are very close to unity. The experimental evidence that the ^^^Ac decay 
chain shows little or no half-life change on cooling has many aspects of which three are 
considered. The 1 K and ~10 mK spectra in Fig. [11] were taken for equal times. There is 
clearly no overall change in activity on cooling from 1 K to the lowest temperature, although 
anisotropy in the ^^^Bi transitions is readily seen. 

(i) The close similarity of the alpha spectra at room temperature, 1 K and 10 mK, Figs.fTOlfTT] 
shows that the predictions of the Full scenario introducing new alpha transitions and elim- 
inating others, did not occur. The ^^^Th, ^^^Ra and ^^^Rn peaks were predicted to vanish 
before reaching 1 K. 

(ii) The observed lack of anisotropy in decay of ^^^Th (spin 1/2), ^^^Ra (weak hyperfine 
interaction) and ^^^Po (extremely short half-life compared to relaxation time) are all in 
agreement with predictions of theory for the reasons cited [13]. These results would be 
strongly distorted in the Partial scenario by the 30% increase in activity on cooling from 1 
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FIG. 10: Room temperature spectrum of ^^^Ac chain alpha decay. 

K to 25 mK. 

(iii) The observed large anisotropies of the 6623 keV alpha transition from ^^^Bi (see Fig. [T2|) 
have small experimental errors. They show the correct relationship in two directions, 86° 
and 24°, one requiring increase and the other decrease of the recorded activity, and strong 
temperature dependence leading to constant maximum values at the lowest temperatures 
reached. As discussed above for the ^^^Bi 292 keV transition in the ^^^Ra chain, all these 
aspects of the data would be disturbed by any increase in activity between 1 K and the 
lowest temperatures caused by changing lifetime. For this case the precision of each data 
point in Fig. [T2]is The data were fitted to give the magnetic dipole moment of ^^^Bi 

dingthe value | yU | = 3.79(7) n.m. in excellent agreement with theoretical prediction 



yie. 



19|, 



20l |. The predicted upward shift of all data in this figure in the Partial scenario is 
close to 30% and would be impossible to miss. Thus the limit on half-life change from the 
mass chain is 1% or lower, based on the lack of activity increase in cooling from IK 
to ~25 mK, on the constancy of (normally corrected) warm counts and detailed evaluation 
of data on many transitions. The absence of change in the alpha spectrum, although well 
established qualitatively, places a less stringent limit on possible half-life change since to be 
clearly observed it requires considerable increase in the alpha decay branching from ^^^Ac. 
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FIG. 11: ^^^Ac decay chain: Alpha spectra in two detectors, upper panels at 86", lower at 24", to 
the polarization axis. The right-hand panels show spectra taken at 1 K, left-hand panels at ~10 
niK. 
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FIG. 12: ^^'^Ac decay chain: ^^^Bi 6623 keV alpha transition anisotropy vs inverse temperature J] 
VII. DISCUSSION AND CONCLUSIONS 



In this paper evidence has been presented concerning the prediction that introducing 
alpha activities into metals and cooling them to temperatures of a few Kelvin and below 



23 



should produce major changes in their decay rates and the appearance of their alpha spectra. 
On this evidence, no such changes have been observed with an upper level of one percent, 
at 1 K, a factor of about 10^ smaller than the predicted changes. Two of the three nuclear 
moments deduced from these experiments, which depended upon correct source decay treat- 
ment, have been confirmed by independent subsequent measurements; the third, measured 
solely by this technique, agrees well with theoretical calculation. 
The evidence presented has been of five types, namely 

1. consistency of warm counts corrected using 'normal' half-lives - a 1 K property. 

2. absence of activity change between 1 K and millikelvin temperatures as shown by 
smooth temperature variation of observed anisotropics and the correct low temper- 
ature saturation behaviour of strong anisotropics. Large effects are predicted, even 
in the Partial scenario. However, these predictons depend upon accepting the T~^/^ 
dependence of the screening energy Ud as applicable down to millikelvin temperature 
- the original suggestions in Ref. ij made no limit of temperature range and took as 
their example cooling to 4.2 K. 

3. absence of anisotropy where isotope properties lead to no polarisation. 

4. no loss of peaks or appearance of additional peaks in alpha spectra between room 
temperature and millikelvin temperatures. 

5. evidence that ratios of relaxation time to half-life are not changed enough to prevent 
thermal equilibrium in some shorter lived ground states. 

Types 3-5 are qualitative rather than quantitative features which, whilst denying the pos- 
sibility of the Full scenario are more difficult to use to place specific limits on much smaller 
effects. 

The emphasis in this paper has been on predictions for alpha decay, not only because 
the effects suggested are so large, but also because, if confirmed they would have impor- 
tant consequences with regard to long-lived waste disposal. However the same ideas also 



predict more modest changes in beta decay half lives [2ll]. The LTNO technique has, for 
more than forty years, studied radioactive samples cooled to sub-Kelvin temperatures. Since 
the introduction of methods of cooling of metal samples to millikelvin temperatures in the 
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early 1960 's, with activity melted, diffused and, more recently, implanted into metals, it has 
been a routine element of such experiments to correct for decay in the source activity. Such 
corrections have always used the 'normal' half-life of the isotopes concerned. Many tens of 
different beta-emitting isotopes have been oriented, with lifetimes ranging from many years 
to a few hours and even less. The host metals used in different studies include the more 
common iron, cobalt and nickel but also copper, silver, gold and many lanthanides. No 
evidence for any deviation from the 'normal' lifetime for these isotopes has been reported 
by the many research groups involved in this work (for comprehensive listing see 10| and 
references therein or search for results on low temperature nuclear orientation). Although 
not directly addressed in this paper, a firm upper limit on any lifetime changes associated 
with cooling beta activities in pure metals can be set at or below 1%. (The well estab- 
lished phenomenon of change in electron capture half-lives with variation of inner electron 
concentration in different environments is excluded from this remark.) 



Since Ref. ^ appeared there have been severa 



follow-up papers. Amongst these the 



22 1 claim to have observed a 6% lifetime 



original proponents of the effect, Raiola et al. 
decrease on cooling a sample of ^^°Po, implanted into Cu, to 12K - a result 'consistent in 



22j (i.e. 6% compared to 1000%, or 1 part in 



sign but significantly smaller than expected' 

150 of the prediction). Jeppesen et al. ^ report changes in half-life at room temperature 
between decay of ^^^Fr implanted into metals and non-metals of order 0.4%, with comparable 
errors, where the Full prediction would be reduction by a factor ~2 (1 part in 500 of the 
prediction) . 

The analysis presented in this paper places much more stringent limits on the size of any 
effect. The upper limit of 1% on any effect at 1 K, compared to predictions of factors above 
10^ represents 1 part in 10^. The evidence for lack of effect down to 25 mK would, if the 
proposed theory were simply accepted at all temperatures, place limits lower than 1 part in 
10^0. 

One may conclude that the suggestion of alleviating problems associated with the dis- 
posal of long-lived radioactive waste by cooling in pure metals is unrealistic. Concerns 
expressed in Ref. l|] regarding the reliability of other less topical but scientifically equally 
important situations where reliance is placed upon the use of 'normal' half-lives in a metallic 
environment, are also misplaced. 
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